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I.  IffTRODUCTIO* 


The  work  described  in  this  report  was  motivated  by  the  need  to  explain 
observed  inconsistencies  between  measured  dynamic  response  data  and  response 
predictions  of  various  analytical  and  finite  element  models  of  certain  rod¬ 
like  structures. 

It  was  felt  that  a  sufficiently  general  mathematical  model  of  this  class 
of  structures  might  illuminate  hitherto  unnoticed  dynamic  effects  and 
interactions  caused  by  such  phenomena  as  large  displacements  and  rotations, 
initial  space  curvature,  variable  shape  and  area  of  cross  section,  initial 
twist  of  cross  section  principal  coordinates  along  the  rod  axis,  transverse 
shear  deformation  and  rotatory  inertia. 

Although  the  study  of  vibrations  of  rods  and  bars  is  one  of  the  most 
ancient  in  structural  mechanics,  a  completely  general  formulation  of  the 
problem  is  difficult  to  find  in  the  literature. 

One  of  the  seminal  treatments  of  the  derivation  of  equations  of  motion  of 
rods  with  initial  curvature  appears  in  the  text  on  Elasticity  by  A.E.H. 

Love.  Although  restrictive  forms  of  Love's  equations  are  used  directly,  or 
are  re-derived,  by  subsequent  investigators  of  small  amplitude  vibrations  of 
curved  rods,  there  are  deficiencies  in  Love's  work  which  render  his  strain- 
displacement  equations  unsuitable  as  the  starting  point  for  a  general 
examination  of  the  motion  of  rods  with  space  curvature. 

In  Love's  work,  the  final  state  of  deformation  is  assumed  to  be  such  that 
cross  sections  remain  plane  and  normal  to  the  centerline  of  the  deformed 
rod.  By  these  assumptions  not  only  is  transverse  shear  deformation  excluded 
but,  as  will  be  shown,  kinematical  inconsistencies  are  introduced. 

Love  employs  a  two  stage  deformation  process  consisting  of  first  the 
imposition  of  twist  and  curvature  "n  the  initially  straight  rod  centerline 
followed  by  imposition  of  small  displacements  to  points  of  the  cross 
section.  This  approach  leads  to  confusion  of  local  coordinate  systems  and  of 
displacement  components  which  renders  suspect  his  expression  for  the  change  in 
length  of  line  segment  connecting  adjacent  points.  The  generality  of  the 
strain-displacement  equations  presented  is  further  limited  by  the  assumptions 
that  displacements  are  small  and  that  dimensions  of  the  rod  cross  section  are 
small  in  comparison  to  the  radius  of  curvature.  Finally,  strain-displacement 
equations  for  an  initially  curved  rod  appear  explicitly  only  in  a  restrictive, 
plane  curvature,  form. 

Equations  governing  both  static  and  dynamic  deformation  of  curved  rods 
are  employed  in  such  areas  of  structural  mechanics  as  vibrations  of  curved 
beams,  large  deflection  and  stability  analysis  of  beams  and  columns  and  non¬ 
linear  dynamics,  including  stability  of  oscillations  of  beams  and  strings.  A 
bibliography  of  this  literature  with  emphasis  of  curved  beam  vibrations  is 
included  with  this  report. 


A.E.H.  Love,  A  Treatise  On  The  Mathematical  Theory  of  Elasticity ,  4th  Bi., 
Dover  Publication,  New  York,  192?. 


In  spite  of  the  large  number  of  papers  dealing  with  these  subjects,  there 
does  not  appear  to  have  been  a  re-examination  of  Love's  strain-displacement 
equations  or  a  more  general  independent  derivation.  The  equations  of  motion 
employed  by  the  various  investigators  are  either  those  presented  by  Love  or 
are  derived  for  each  particular  application  on  an  "ad-hoc"  basis. 

In  this  report  a  derivation  of  the  equilibrium  equations  for  a  rod  with 
general  space  curvature  is  presented  in  Section  II.  This  derivation  is 
similar  to  that  presented  by  Love,  or  more  recently  by  Kachanov,  although  it 
i9  somewhat  more  rigorous  than  those  in  its  examination  of  the  relative 
importance  of  higher  order  terms  in  Taylor's  series  expansions. 

A  completely  general  set  of  small  strain  but  large  displacement  and 
rotation  strain-displacement  equations  is  then  formulated  in  Section  III  for  a 
rod  with  initial  twist  and  arbitrary  space  curvature.  In  these  equations  the 
displacement  components  are  referred  to  the  local  tangent,  principal  normal 
and  bi-normal  space  curve  coordinates  rather  than  to  fixed  global 
coordinates.  No  assumptions  regarding  the  structural  action  or  mode  of 
deformation  of  the  rod  are  made  so  that  the  functional  form  of  the 
displacement  components  is  left  unspecified. 

In  Section  IV  a  linearized  form  of  the  general  strain-displacement 
equations  is  first  presented.  A  technical  theory  of  rods  with  space  curvature 
is  then  developed  based  in  the  usual  Learn  theory  assumptions  that  cross 
sections  remain  plane  and  undeformed.  The  displacements  are  then  expressed  as 
appropriate  linear  combinations  of  three  central  curve  displacement  components 
and  three  rotation  components. 

Based  upon  these  displacement  functions,  the  inertia  force  and  moment 
terms  appearing  in  the  equilibrium  equations  can  then  be  formulated  in  terms 
of  time  derivatives  of  these  six  displacement  variables. 

The  linearized  strain  displacement  equations  phrased  in  terms  of  the 
central  curve  displacements  and  cross  section  rotations  are  then  used  to 
derive  force-moment-strain-curvature  equations  for  a  linear  elastic  rod.  The 
equations  are  then  further  simplified  by  elimination  of  transverse  shear  and 
rotatory  inertia  effects.  When  these  equations  are  combined  the  equilibrium 
equations,  a  system  of  four  coupled  equilibrium  equations  in  the  four 
displacement  and  rotation  variables  results. 

The  final  section  of  the  report  is  devoted  to  examination  of  the 
displacement  equilibrium  equations  for  special  cases  in  order  to  examine  the 
conditions  under  which  various  uncouplings  of  motions  occur. 


L.M.  Kachanov,  "A  Brief  Course  On  The  Theory  of  Buckling’’  t&AM  Report  No. 
441.  University  of  Illinois  Urbana-Champaign ,  1980. 

8 


II.  EQUILIBRIUM  EQUATIONS  FOR  A  BOD  ELEMENT 


A.  Force  gad  Maggot  Equilibrium 

A  typical  element  of  a  rod  with  space  curvature  is  shown  in  Figure  1. 

The  line  joining  centroids  of  cross  sections  is  a  space  curve  whose 
radius  (p),  cr  curvature  tc*l/p,  and  torsion  (X)  are  arbitrary  functions  of 
arc  length  (s).  The  shape  of  the  cross  sections  is  arbitrary  but  assumed  to 
be  a  continuous  function  of  s.  At  thf  centroid  of  any  cross  section  an 
orthogonal  coordinate  system  (x,y,z)  is  constructed  with  corresponding  unit 
vectors  i,  j,  k,  such  that  k  is  the  unit  tangent  vector  to  the  space  curve 
determined  by  the  line  of  centroids.  Vector  k  therefore  points  in  the 
direction  of  increasing  s.  Unit  vectors  i  and  j  are  the  principal  normal  and 
bi-normal  vectors  respectively  so  that  jxk  *  i. 

The  principal  directions  of  the  cross  section  area  of  the  rod  with 
respect  to  the  centroid  are,  in  general,  inclined  with  respect  to  the  x,  y 
coordinates.  The  angle  between  these  two  sets  of  coordinate  axes  may  vary 
continuously  with  8  for  the  case  of  a  pre-twisted  rod. 

Figure  2  shows  a  section  of  curved  rod  of  incremental  length  6s  lying 
between  points  0  and  O'  of  the  space  curve  of  centroids.  By  the  usual 
definition,  the  face  at  0  is  a  negative  face  since  its  outward  normal  points 
in  the  direction  opposite  to  the  local  tangent  vector  while  the  face  at  O'  is 
a  positive  face.  The  usual  sign  convention  for  stress  components  in  positive 
and  negative  faces  will  be  employed. 


The 
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Appendix  A  presents  a  discussion  of  space  curve  geometry 
Equations  are  numbered  consecutively  within  each  section. 
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Figure  5.  Adjacent  Cross  Sections  Before  Defoliation 
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III.  STSAIH-DISPLACKMEHT  EQOATIOHS 

Equations  relating  the  components  of  strain  at  a  point  in  the  rod  to  the 
displacement  field  are  derived  by  considering  the  change  in  length  of  a  line 
segment  connecting  to  arbitrarily  close  points  of  the  rod  as  the  rod  is 
deformed. 

Figure  5  shows  two  adjacent  cross-sections  of  the  rod  before 

deformation.  The  centroidal  curve  has  initial  curvature  k  and  torsion  A  . 

o  o 

The  centroid  of  one  cross-section  is  point  0  and  the  other,  at  a  small 
distance  6s  along  the  curve,  is  O'.  The  vector  principal  normal,  binormal  and 
tangent  at  0  are  i,  j  and  k  respectively  while  the  corresponding  vectors 
at  O'  are  denoted  by  i,  j  and  k.  Coordinates  x,  y,  z  directed  along  i,  j  and 
k  respectively  comprise  a  local  coordinate  system  with  origin  at  0. 
Corresponding  coordinates  at  O' are  x,  y  and  z. 

Point  Q  lies  in  the  cross-section  at  0  while  Q'  lies  in  that  at  O'.  It 
is  assumed  that  Q  and  Q'  are  initially  arbitrarily  close  so  that  the  position 
vector  of  Q',  r  can  be  expressed  as 


r  *  6x  i  +  6y  j  +  6z  k 


(1) 


where  6x,  6y,  and  6z  are  arbitrarily  small  distances. 

If  n  is  a  unit  vector  parallel  to  r  and  directed  from  Q  to  Q',  and  r  is 
the  magnitude  of  r  .  Then 

r*marili  +  rmj  +  rnk  (2) 

when  Jl,m,  and  n  are  the  direction  cosines  of  n  relation  to  i,  j  and  k 
respectively.  Equations  (1)  and  (2)  imply  the  relationships 

6x  *  r£,  6y  a  rra,  and  6z  =  rn  .  (3) 

The  length  of  the  position  vector  joining  material  points  initially  at  Q 
and  Q'  as  the  rod  undergoes  a  deformation  which  carries  these  points 
to  Q  and  Q'  respectively  is  next  determined. 
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For  a  discussion  of  various  representations  of  the  angular  acceleration  and 
angular  velocity  vectors  the  reader  is  referred  to  a  recent  book  by  Kane,  et. 
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Upon  integrating  over  the  area  of  the  cross  section,  A,  the  inertia  forces  per 
unit  length  of  rod  arising  from  the  three  body  force  components  are  next 
obtained. 
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Finally  the  inertia  force  couples  are  presented  in  terms  of  the  three 
acceleration  components 
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T.R.  Kane,  P.W.  Likins ,  and  D.A.  Levinson,  " Spacecraft  Dynamics,"  MoGraw 
Hill ,  New  York,  1983. 


The  reader  is  referred  to  a  paper  by  E.  Reissner  for  corresponding  large 
displacement  curved  beam  equilibrium  equations  in  which  the  force  and  moment 
components  are  referred  to  the  undeformed  coordinates. 

B.  Inertia  Forces  and  Couples 

The  inertia  body  force,  Fj,  acting  at  any  point,  Q,  in  a  cross  section  is 
defined  by 


~  Y 


where  y  is  the  mass  density  of  the  material  of  that  point,  uQ  is  the 
displacement  vector  associated  with  that  point  and  superscript  "dot"  indicates 
differentiation  with  respect  to  time. 

The  displacement  vector  Uq  is  expressed  in  terms  of  the  local  space  curve 
coordinates  as 


ui  +  vj+wk 


where  u,  v,  and  w  are  the  scalar  displacement  components  in  the  x,  y  and  z 
directions  respectively. 

'Since  the  local  coordinate  system  is  moving  in  the  fixed  (inertial) 
reference  frame,  the  acceleration  vector  may  be  expressed  in  the  form: 


Uq  -  a^  +  o  x  Uq+  mx(tdxuq)  +  2  u  x 


where 


ui  +  vj+wk, 


^*ui  +  vj+wk, 

and  m  and  a  are  the  angular  velocity  and  angular  acceleration  vectors 
respectively  of  the  local  coordinates  in  the  fixed  reference  frame.  For 
problems  involving  small  rotations  and/or  small  displacements  of  the  rod  cross 
the  non-linear  terms  in  Eq.  (23)  may  often  be  neglected  in  comparison  with  «q. 

Since  this  report  will  be  concerned  primarily  with  equations  governing 
small  motions  the  non-linear  terms  in  Eq.  (23)  will  henceforth  be  neglected. 


E.  Reissner,  "On  One-Dimensional  Large-Diaplacement  Finite-Strain  Beam 
Theory,"  Studies  in  Applied  Math »  Vol  52,  pp  87-95 ,  1973. 


Although  this  manner  of  derivation  of  the  equations  of  force  and  moment 
equilibrium  for  a  curved  rod  used  above  is  different  from  that  presented  by 
Love1  the  resulting  equations  are  essentially  identical.  Differences  in  form 
result  from  differences  in  force  and  moment  sign  conventions  and  because  Love 
employs  local  coordinate  axes  which  are  the  principal  directions  of  the  area 
of  the  local  cross  section  of  the  rod  rather  than  the  natural  space  curve 
coordinates. 

In  general,  the  local  x,  y,  z,  coordinates  are  the  deformed  coordinates 
of  the  centroidal  curve  of  the  rod.  The  torsion  and  curvature  terras  appearing 
in  Eqs .  (12)  and  (21)  are  the  total  values  of  these  quantities  consisting  of 
the  suras  of  their  initial  values  and  their  respective  changes  due  to 
deformation  of  the  centroidal  curve.  For  small  deformation  analysis  the 
initial  coordinates  and  initial  values  of  \  and  k  are  employed. 
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+x  .  av. 


£'  x  V  (a  +  6a)  -  [k  6a  +  k/2  i  (6a)‘]  x  [▼  (a  )  +  ff)  |  +6 


+  77  I  +  V2«s>2] 

08  8 


k  x  V  (s+)  6s  +  [V2k  i  x  V(s  +  )  +  k  x  -|j  |  +]  (6s)2  and  (17) 

8 


2 

M(s  +  6s)  -  M(s+)  +  |f  |  +«*  +  1/2  |  +(6a)2  ,  (18) 

8  3s  8 


where  again  V(s+)  and  M(s+)  indicate  the  force  and  moment  resultants  acting 
on  the  positive  face  at  6s  *  0  so  that  M(s+)  *  -  H(s). 

Finally,  substitution  of  the  result  given  by  Eqs.  (14)  and  (16)  through 
(18)  into  Eq.  (13)  and  division  by  6s  yields  the  following  vector  moment 
equilibrium  equations  in  the  limit  as  6a+0. 


+  +  kxV(s+)  +  *(s) 


To  obtain  the  scalar  equations  represented  by  Eq.  (19),  M(s  )  is  expressed  in 
terms  of  its  vector  components  and  the  Serret-Frenet  formulae  are  used  to 
transform  the  derivatives  of  the  unit  vectors 


3M  3M  3M 

j-*t  *  *  Jk  ♦  k  ii  .  nil  *  1  ? 

3s  3s  os  x  ds  y  ds  z  ds 


3M  3M  3M 

i  +  r-2-  j  +  t-5-  k  +  M  (Xj  -  Kk)  +  M  (-A)i 

38  38  3s  x  ■*  y 


+  <  M  i  . 

z 


Also,  since 
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In  Figure  4,  R0,  R'0  and  P  are  position  vectors  of  points,  at  positions 
"s",  "s  +  6s”  and  "s  ♦  C"  relative  to  a  fixed  point  Q.  If  R'(s)  is  the 
position  vector  of  any  point  in  the  space  curve  relative  to  0  we  nay  write 


R’  -  R 

o  o 


dR/ds  I  6s  +  4-  ^-y  I  (6s)2  + 

,s  2  ds2  8 


2 

dR/ds  I  6s  +  4-  ^-y  I  (6s)2  +  ... 

1 8  2  .  Z  1 8 

ds 


k6s  +  k/2  i  (6s)  , 


(15a) 


since  by  definition  k  “  dR/ds  and  by  Eq.  (10a)  dk/ds  "  k  i. 
In  a  similar  way  it  is  found  that 


P  -  R 


kC  +  k/2  i  c' 


(15b) 


Combining  Eqs.  (7)  and  (15b)  yields 


R  x  T  -  k(s)  x  T(s)  5  +  [</2  ixT  +  kx|j|^2  + 


Therefore. 


/  RxTdj  *  kxl|  id?  +  [*c/2  ixT  +  k  x  3T/3s]  /  Oc 


-  k(s)  x  T  V2  (6s)2  . 


Next, 
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av 

_ X 

38 


X  V  )i 

y 


0 


(11) 


av 

(? r  -  KVx^ k  +  (tev  z  +  - 


+  (X  V  +  3V  /3s) j  +  Ti  +  Tj  +  Tk- 
x  y  x  yJ  z 

Finally,  Eq.  (11)  may  be  written  as  three  scalar  equations  which  express  the 
condition  of  force  equilibrium  for  a  beam  element 


a  v^/a s 

a  v  /as 
y 

a  v  /as 

z 


-XV  +  kV  +  T  -  o  , 
y  z  x 

+  XV  +  T  ■  0  ,  and 

x  y 


kV  +  T  -  0  . 

X  z 


(12a) 

(12b) 

(12c) 


Next,  the  conditions  for  moment  equilibrium  are  formulated  from 
consideration  of  the  free-body  diagram  of  Figure  4,  Summing  moments  about 
point  0  yields 


M(s)  +  M(s  +  6s)  +  R'xVCa  +  6s) 

8+&8  S+6s 

+  /  *xT(s)ds  +  /  ■  ds  -  0  (13) 

8  8 


Again  the  terms  appearing  in  Eq.  (13)  must  be  expressed  in  terms  of  their 
values  and  those  of  their  derivatives  at  'V 


m(s  +  ?) 

s+6  s 

/  *ds  * 


,<*)  *  If  I.  5  *721. 

3s 


2 

C  /2  +  ...  and 

s 


6s 


6s 


n(s)  J  dc  +  |8  /  4dC  ♦ 

o  o 


1  I  f63  2 

2  .  2  1 8  C  dC  +  **• 

3s  o 


"  "(8)6s  +  If  I 

s 


(6s)2  + 

2 
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(14) 


8+68 

/  Tdi 
8 


T(a)6s  +  (6s)2/2  ♦  ... 

08  8 


(8) 


Returning  to  Eq.  (7),  we  next  express  V(a  +  6s)  in  terns  of  its  value 
and  those  of  its  derivatives  at  is  ■  0  using  vector  Taylor's  series. 


?(s  +  68)  ■  V(s+)  +  |  +  6s  +  ...  (9) 

s 

The  symbol  "s+"  is  used  in  Eq.  (9)  as  a  reminder  that  this  force  acts  on  a 
positive  face  of  the  element.  Therefore, 


V(s+ ) 


V  i  + 
x 


V  j 

y 


+  V  k  and 
z 


(9a) 


3V 

9  s 


+ 


V  ~ 
x  ds 


+  V 

y 


V 

z 


dk 

ds  * 


(9b) 


j 


Expressions  for  the  rate  of  change  of  the  unit  tangent,  principal  normal 
and  binormal  vectors  with  respect  to  arc  length  are  given  by  the  Serret  - 
Frenet  formulae  (Appendix  A)  as 


dk/ds 

ii  . 

ds 

di/da 


1/pis  k  i  , 


X i  ,  and 


Xj  -  xk 


(10a) 

(10b) 


(10c) 


where  k(*  1/p)  is  the  local  curvature  of  the  line  of  centroids.  Upon 
substitution  of  Eq.  (8)  through  Eq.  (10)  into  (6)  and  division  by  6s,  the 
following  vector  equation  results  in  the  limit  as  6s  approaches  zero, 


•XL, 
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The  surface  tractions  and  body  forces  may  also  contribute  to  a  moment 
resultant  per  unit  length  about  the  three  coordinate  axes.  This  moment 
resultant  is  defined  by 


m  k  . 
z 


For  problems  involving  motion  of  the  rod  element,  the  force  T  will 
contain  terms  involving  inertia  forces  while  the  moment  ■  may  contain  moments 
of  inertia  forces  about  the  three  axes. 

A  free-body  diagram  of  a  rod  element  is  shown  in  Figure  4.  In  Figure  4, 
points  0  and  O'  are  again  points  in  the  curve  of  centroids  at  positions  "s" 
and  "s  +  6s"  respectively.  The  force  and  moment  resultants  acting  on  the  rod 
cross  sections  of  these  positions  are  indicated  as  well  as  the  surface 
traction  vector  and  surface  traction  moment  resultant  at  a  typical  point 
s  +  C  where  s<s  +  C<s  +  6s. 

_  t  I 

R  is  the  positiont vector  of  the  point  x  +  C  relative  to  0  while  R  is 
the  position  vector  of  0  relative  to  0. 

The  condition  of  force  equilibrium  applied  to  the  rod  element  yields  the 
following  equation 


V(s)  +  V(s  +  6s) 


s+6s 
/  T  ds 


1  32T 


«■♦«>  -  «•>  ♦  Sf  l.<  ♦  i  H  l.<‘ 

os 


where 


<(s  +  O  <(s  +  6s)  and  ds  *  d£  . 


Therefore, 


/  T  ds  -  T(s)  /  de  ♦  P  J  4 d; 


9 

I  O 

+  o  o 


/  <Td c  +  ... 


(2c) 


M  ■  f(x  a  -  y  o  )dA  . 
z  i  yz  ’  xz 

Based  on  the  above  definitions,  the  vector  force  and  moment  resultant  of 
the  tractions  acting  on  the  face  with  centroid  "0"  at  point  "s"  on  the 
space  curve  of  centroids  is 


V(s)  ■  -  V  (s)i(s)  -  V  (s)j(s)  -  V  (s)k(s)  and  (3a) 

x  y  x 


M(s )  =  -  M  (s)i(s)  -  M  (s)j(s)  -  M  (s)k(s) 

x  y  z 


(3b) 


Similarily,  these  vector  quantities  acting  at  s  +  6s  on  the  face  with  centroid 
O'  are 


V(s  +  6s) 


M(s  +  6s) 


=  V  (s+6s)i(s  +  6s)  +  V  (s  +  6s)j(s  +  6s) 

x  y 

+  V  (s  +  6s)k(s  +  6s)  and 
z 

=  M  (s  +  6s)i(s  +  6s)  +  M  (s  ♦  6s)j(s  +  6s) 

x  y 

+  M  (s  +  6s)k(s  +  6s)  . 
z 


(4a) 


(4b) 


These  quantities  are  illustrated  in  Figure  3. 

Note  that  the  unit  vectors  used  to  describe  the  directions  of  the  vector 
forces  and  moments  at  s  are  different  from  those  at  s  +  6s  since  the 
directions  of  the  space  curve  unit  vectors  are  constantly  changing  along  the 
curve . 


The  resultant  per  unit  length  of  the  tractions  acting  on  the  lateral 
surface  and  of  the  body  forces  acting  in  the  interior  of  the  element  is  given 
by 


T*Ti  +  Tj  +  Tk  (5a) 

x  yJ  z 

where  T  ,  T  ,  T  have  the  dimensions  of  force  per  unit  length  and  i,  j,  k  are 
x  y  z 

local  unit  vectors. 


Figure  6  illustrates  the  initial  and  final  positions  of  these  material 
points  as  Q  and  Q*  undergo  displacements  and  u-,.  The  magnitude  of  the 
position  vector  of  Q'  with  respect  to  Q,  r  is  found  by  use  of  the  following 
vector  equation  relating  the  initial  and  final  relative  position  vectors  of 
these  points  to  their  respective  displacement  vectors. 


r  *  r  + 


V  -  "Q  • 


(4) 


The  displacement  components  u,  v,  and  w  at  point  Q  are  again  given  by 


ui  +  vj  +  wk. 


(5) 


The  displacement  vector  at  point  Q',  adjacent  to  point  Q,  can  then  be 
expressed  as 


u, 


V 


*  ? 


6x 


3y 


6y 


6z 


(6) 


where  higher  order  terms  in  this  Taylor's  series  expansion  of  the  displacement 
vector  about  point  Q  have  been  dropped.  The  vector  derivatives  appearing  in 
Eq.  (6)  are  next  presented  in  terms  of  their  vector  components. 


3u, 


l!a 

3y 


3z 


3v 


3w 


Q  8u  ........  . 

3lT  3x  1  3x  J  3x  k  ’ 


(7a) 


3u  .  .  3v  .  .  3w  . 

37  1  *  37  1  *  37  k  > 


(7b) 


3u  .  .  3v  .  3w  .  3  i  3j  3k 

371  +  TTJ  +  37k  +  u3l+v3;+w37 


.  .  3v  .  A  3w 

i  1  j  +  r 


3u 


k  +  u  (X  j  -  x  k) 


vX  i  +  w  k 


■j— ^  “  (t~  -  Xv+Kw)i  +  (~  +  X  u)  j  +  (|^  -  utc  )  k 

3z  '“3*  o  o'  '“Sz  o  '  "*  v3z  o' 


where  the  Serret-Frenet  formulae  have  been  used  to  express  the  rate  of  change 
of  the  local  coordinate  unit  vectors  along  the  curve  in  terms  of  the  local 
curvature  and  torsion. 

Substitution  of  Eqs.  (5)  through  (7)  into  (4)  and  using  Eqs.  (3)  yields 
the  following  result  for  the  position  vector  of  Q1  relative  to  Q: 


r./,  .  3uv  .  3u  ^  (3u 

r  U(l  +  d  +  ®  r  +  n  b 


X  V  +  K  w)]  i 
o  o 


+  rU  If  +  “U+fj)  +  +  3  i 


+  r[Jl  +  m  -r^-  +  n(l-<u  +  ]  k  . 

3x  3y  v  o  ozJ 


If  the  magnitudes  of  r  and  r  a’-e  denoted  by  r  and  r  respectively,  then  the 
unit  elongation,  e,  of  the  original  line  segment  is  defined  by 


r  -  r(l  +  e)  . 


If  e  is  a  small  quantity,  then 


r  2  -  r2(l  +  2e) 


(10a) 


v2  (4  - 1)  • 


(10b) 


Upon  squaring  the  magnitudes  of  the  vector  r  given  by  Eq.  (8), 
substituting  the  result  into  Eq.  (10b),  and  grouping  the  resulting  terms  as 
coefficients  of  the  products  of  the  direction  cosines  Z ,  m  and  n,  the 
elongation,  e,  may  be  expressed  in  the  following  form, 
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By  definition,  the  quantities  e  ,  e  ,  e  ,  e  ,  e  .  and  e  are 

xx  yy  zz  xy  .xz  yz 

and  shear  strains  at  Q  referred  to  the  local  coordinates  at  ?.  Gqs 
therefore  constitute  the  required  strain-displacement  equations. 


(11) 

(12) 

(12b) 

(12c/ 

(12d) 

(12e) 

( 1 2  f  ) 

the  normal 

.  (12) 
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Clearly,  various  approximate  forms  of  Eq.  (12)  may  be  formulated 
depending  upon  the  simplifying  assumptions  appropriate  for  a  given  problem. 
These  formulations  may  postulate  small  deformations  with  large  displacements; 
small  deformation  with  small  displacements  but  large  curvatures  and/or 
torsions;  small  curvatures  and  torsions;  or  various  other  combinations  of 
simplifying  assumptions.  The  reader  may  verify  that  the  technical  theory  of 
curved  beams,  in  which  the  normal  strain  is  a  non-linear  function  of  the 
transverse  coordinate,  requires  retention  of  terms  involving  squares  of 
curvatures  and  of  displacements.  It  is  also  noted  that  consideration  of 
buckling  or  of  stretch-stiffening  effects  requires  retention  of  non-linear 
terms  in  Eq.  (12). 

The  simplest  form  of  strain  displacement  equations  which  contains  the 
effect  of  initial  curvature  and  torsion  of  the  centroidal  curve  is  obtained  by 
dropping  all  non-linear  terms  in  displacements  and  displacement  gradients. 

The  strain-displacement  equations  then  become 


xx  dx  » 


(13a) 


(13b) 


V  * 


(13c) 


Exy  "  *  *  £)  • 


(13d) 


e  ■  V2  +  3—  -  A  v  +  k  w)  ,  and 

xz  z  '  3  z  3x  o  o' 


(13e) 


(13f) 


At  this  point  the  functional  form  of  the  displacement  components  in  terms  of 
the  coordinates  is  completely  unspecified.  In  the  next  section  the  strain- 
displacement  equations  are  put  in  a  form  which  corresponds  to  the  usual 
technical  beam  theory. 

IV.  A  TECHNICAL  THEORY  OF  RODS  WITH  INITIAL  CURVATURE  AND  TORSION 
A.  Displacement  Functions 


A  linearized  "engineering"  theory  of  curved  rods  is  next  developed  on  the 
basis  of  assumptions  employed  in  the  strength  of  material  formulations  for 
the  bending  and  torsion  of  prismatic  bars. 
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The  fundamental  assumption  of  both  technical  beam  theory  and  that  of 
torsion  of  straight  rods  is  cross  sections  unchanged  in  size  and  shape  after 
the  structure  is  bent  or  twisted.  It  is  not  necessary  to  initially  postulate 
that  cross  sections  remain  normal  to  the  bent  centerline. 


A  typical  cross  section  is  therefore  assumed  to  undergo  a  small  rotation 
about  each  of  the  three  coordinate  axes  while  its  centroid  is  displaced  along 
each  of  these  axes.  These  assumptions  imply  the  following  form  for  the 
displacement  functions, 


w(x,y,z) 

v(x,y,z) 

u(x,y,z) 


x  a  (z)  +  y  a  (z)  , 

(la) 

y  J  x 

x  ( z )  ,  and 

(lb) 

y  4>(z)  . 

(lc) 

In  Eq.  (1),  uQ,  vQ  and  wQ  are  the  displacement  components  at  a  point  on  the 
centroidal  curve  (the  origin  of  the  local  x,y  coordinates)  while 
a  ,  a  and  <j>  are  small  angles  of  rotation  about  the  x,  y  and  z  axes  in  that 
or  <fer  ,yaa  shown  in  Figure  7. 

It  is  noted  that  theory  of  torsion  of  bars  with  non-circular  cross 
sections  requires  the  addition  to  Eq.  (1)  of  a  term  representing  warping  of 
the  cross  section.  Introduction  of  such  a  term  would  greatly  complicate  the 
ensuing  development  with  very  little  effect  on  the  state  of  stress  other  than 
near  a  fixed  boundary.  For  this  reason,  the  displacement  functions  of  Eq.  (1) 
will  be  employed  irrespective  of  the  shape  of  the  cross  section  of  the  bar. 

Warping  of  compact  cross  sections  increases  as  the  ratio  of  the  principal 
second  moments  of  area  differs  increasingly  from  unity.  The  warping 
phenomenon  is  thoroughly  treated  in  texts  by  Timoshenko  and  Goodier  and  by 
Sokolnikof f? 

B.  Strain-Displacement  Equations 

The  linearized  form  of  the  strain-displacement  equations  derived  in 
Section  III  are  used  as  the  basis  for  the  strain-displacement  equations  for 
the  technical  theory. 


b  S.  Timoshenko  and  J.N.  Goodier ,  Theory  of  Elasticity ,  2nd  Ed.,  McGraw- 
Hill,  New  York,  1951. 

$  I.S.  Sokolnikoff ,  Mathematical  Theory  of  Elasticity,  2nd  Ed.,  McGraw-Hill , 
New  York,  1956. 
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Substitution  of  Eq.  (1)  into  Eq.  (III-13)  yields  the  following 
expressions  for  the  components  of  strain. 


e  -  e  ■  e  *  0  , 
xx  yy  xy 


(2a) 


3w  3a  3a 

e  -^-x^+y^  -  <Q(u0  -  y*)  , 


'zz  3z  3z 

3u 


(2b) 


2  e  *  ~r — —  -  a  -  X  v  +  <w  +  y(<  a  -  -x(X<i>  +  tca)  ,  and 

xz  3z  y  o  o  o  o  J  o  x  3z'  o  o  y 

(2c) 


3v  a , 
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Eq.  (2a)  represents  the  condition  that  cross  sections  do  not  deform  and 
thereby  confirm  the  satisfaction  of  the  initial  displacement  assumptions. 

Eqs.  (2c)  and  (2d)  show  that  transverse  shear  deformation  cannot  be  made 
to  vanish  for  all  values  of  x  and  y.  Introduction  of  the  angle  <|>  to  represent 
torsional  rotation  would  lead  to  this  result  even  for  a  straight  bar.  For  the 
curved  bar  it  is  noted  that  the  presence  of  initial  torsion  and  curvature 
prevent  these  strains  from  vanishing  everywhere  even  with  the  omission  of  the 
rotation  <b . 

The  condition  that  cross  sections  remain  normal  to  the  centroidal  curve 
can  be  invoked  locally  at  x*0,  y*0  or,  alternatively,  in  an  average  sense  by 
integrating  Eqs.  (2c)  and  (2d)  over  the  cross  section.  In  either  case,  the 
following  equations  expressing  the  vanishing  of  transverse  shear  deformation 
(in  the  above  meaning)  result: 
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Another  simplifying  assumption  often  invoked  is  that  the  centroidal  curve 
(x*y**0)  remains  unchanged  in  length  during  the  motion  of  the  rod.  From  Eq. 
(2b),  this  condition  may  be  expressed  as 


C.  Inertia  Force*  and  Mowenti 


The  inertia  force  resultants  and  couples  per  unit  length  are  obtained  by 
introducing  the  assumed  displacement  functions  into  Eqs.  (11-25)  and  (11-26) 
and  carrying  out  the  indicated  integrations.  It  should  be  remembered  that  the 
origin  of  the  coordinate  system  is  the  area  centroid  of  the  cross  section  so 
that  integrals  of  the  form  /  xdA  and  /  ydA  vanish  .  The  resulting 

A  .  A 

expressions,  assuming  no  variation  of  mass  density  within  a  cross  section,  are 
presented  below. 
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I  and  I  are  the  second  moments  of  area  of  the  cross  section  for  the 
centroid  about  the  x  and  y  axes  respectively;  IXy  is  the  product  of  area  with 
respect  to  the  centroid  for  the  x  and  y  directions;  and  Jp  is  the  polar  moment 
of  area  for  the  z  direction  with  respect  to  the  centroid. 

The  area  and  shape  of  the  cross  section  of  the  rod  may  vary  in  an 
arbitrary  way  along  the  rod.  The  orientation  of  the  principal  directions  of 
the  cross  section  area  with  respect  to  the  plan  of  curvature  also  may  vary. 

The  quantities  defined  by  Eq.  (6)  should  therefore  be  considered  to  be 
functions  of  the  z  coordinate. 

D.  Moment-Curvature  and  Force-Displacement  Equations 


The  stress-strain  relations  which  will  be  used  for  the  derivation  of  the 
equations  relating  the  force  and  moment  resultants  acting  as  a  cross  section 
of  the  rod  on  one  hand,  to  the  displacement  and  rotation  components,  on  the 
other,  are  predicated  on  the  assumptions  that  the  rod  is  constructed  from  a 
linear  elastic  material  which  is  elastically  homogeneous  in  each  cross  section 
and  that  the  normal  stress  components  in  the  directions  transverse  to  the 
centroidal  curve  are  small  in  comparison  with  the  normal  stress  component 


The  applicable  stress-strain  equations  then  become 
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The  required  equations  are  next  obtained  by  combining  Eqs.  (2)  and  (3)  with  the 
definitions  of  the  various  force  and  moment  resultants  given  by  Eq.  (II-l), 
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The  coefficients  kx  and  ky  appearing  in  Eqs.  (8d)  and  (8e)  are  called 
Timoshenko  shear  coefficients.  These  are  dimensionless  quantities,  dependent 
upon  the  shape  of  the  cross  section,  which  are  introduced  to  compensate  for 
the  fact  that  not  only  is  the  shear  stress  non-uniform  over  the  cross  section 
of  the  rod  but  its  actual  distribution  is  inconsistent  with  the  function  forms 
describing  the  shear  strain-displacement  equations.  Values  of  the  Timoshenko 
shear  coefficients  are  documented  for  various  cross  sections  as,  for  example, 
by  Cowper . 

The  complete  system  of  equations  constituting  the  linear,  technical, 
theory  of  rods  with  space  curvature  is  comprised  of  Eqs.  (11-12)  and  11-21) 
along  with  Eqs.  (4),  (5)  and  (8)  of  this  section.  For  convenience  of  the 
following  exposition,  Eqs.  (11-12)  and  (11-21)  combined  with  (4)  and  (5)  are 
presented  below  where,  in  a  manner  consistent  with  the  employment  of 
linearized  strain-displacement  equations,  the  torsion  and  curvature  are  taken 
as  those  of  the  undeformed  rod. 
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G.R.  Cowper,  "The  Shear  Coefficients  in  Timoshenko  Beam  Theory,"  ASME  J. 
Applied  Mechanics,  Vol  33,  Series  E ,  No.  2,  pp  335-340 ,  1966. 
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In  Eq.  (9),  the  terms  Txs,  TyS,  Tzg  represent  the  portion  of  the  surface  and 
body  forces  acting  which  are  not  acceleration  dependent  while  a..,  m  and  m 

XS  Vo  ZB 

are  the  corresponding  parts  of  the  couples  of  these  forces. 

To  determine  the  motion  of  a  rod  with  space  curvature  using  this 
linearized  formulation,  the  system  of  twelve  coupled  equations  consisting  of 
Eqs.  (8)  and  (9)  above  must  be  solved  for  the  six  force  and  moment  resultants 
and  the  six  displacement  and  rotation  components  subject  to  appropriate 
initial  and  boundary  conditions.  This  formulation  retains  the  effects  of 
transverse  shear  deformation,  rotatory  inertia,  and  extension  of  the 
centroidal  curve.  Since  some  or  all  of  these  effects  can  safely  be  neglected 
for  problems  involving  thin  rods  and  low  frequencies,  more  simplified 
formulations  will  next  be  considered. 

It  is  noted  that  no  restrictions  have  yet  been  imposed  on  the  shape  of 
the  cross  section,  the  orientation  of  the  principal  axes  of  the  cross  section 
with  respect  to  the  local  plane  of  curvature,  or  the  change  in  these 
quantities  along  the  rod.  The  cross  section  area,  the  torsion  and  curvature, 
and  the  tensile  and  shear  moduli  may  also  have  arbitrary  variation  with  z. 
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1.  General  Fora 

In  this  section  the  equilibrium  equations  for  the  rod  are  phrased  in  terms 
of  displacement  functions.  The  effects  of  torsion,  curvature  and  orientation 
of  the  principal  axes  of  the  cross  section  areas  as  well  as  the  variation  of 
these  quantities  along  the  rod  are  then  examined. 

In  the  following  formulation  transverse  shear  strains  are  neglected  in 

the  sense  defined  by  Eq.  (3).  Since  rotatory  inertia  generally  has  less 

effect  upon  predictions  of  frequency  and  motion  than  Jjas  transverse  shear 

deformation,  it  is  consistent  to  set  terms  involving  a  and  a  equal  to  zero 

in  Eq.  (9)  although  rotational  inertia  effects  about  tfie  z  ax?s  are 

retained.  Equations  (8d)  and  (8c)  for  Vx  and  V  respectively  become 

inconsistent  because,  although  and  e  are  set  equal  to  zero  on  the 

average,  equilibrium  considerations  require  non-zero  transverse  shear 

forces.  These  forces  are  eliminated  from  the  equilibrium  equations  by  solving 

for  V  and  V  in  Eqs.  (9d)  and  (9e)  and  substituting  these  results  into  the 
y  x 

remaining  four  equations.  The  equilibrium  equations  then  become 
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In  Eq.  (10)  superscript  "prime"  indicates  differentiation  with  respect  to  z. 
Eq.  (9f)  is  repeated  as  Eq.  (lOd). 
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To  obtain  the  required  displacement  form  of  the  equilibrium  equations,  the 
functions  and  o  are  first  eliminated  from  Eq.  (8)  by  use  of  Eq.  (3).  The 
resulting  expressions  for  the  moments  and  Vz  in  terms  of  the  four  displacement 
and  rotation  functions  are  then  substituted  into  Eq.  (10).  The  resulting 
equations  can  be  presented  in  the  following  form: 
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The  fifty-eight  coefficients  appearing  in  Eq.  (11)  are  presented  in 
Appendix  B.  There  it  is  seen  that  these  coefficients  depend  upon  the 
geometric  properties  of  the  space  curve  described  by  the  centroids  of  the 
cross  sections,  the  geometry  of  the  cross  section  and  its  relationship  to  the 
local  vector  principal  normal  and  binormal,  the  elastic  moduli  of  the 
material,  and  the  variations  of  these  quantities  along  the  rod. 

It  is  noted  that  in  the  general  case  described  by  Eq.  (11),  all  motions 
are  fully  coupled  since  all  variables  appear  in  each  equation.  This  means, 
for  instance,  that  any  exciting  force  will  cause  motion  involving  all  four 
displacement  and  rotation  variables.  Also,  solution  of  the  free  vibration 
problem  will  yield  mode  shapes  which  will  each  have  four  corresponding  natural 
frequencies . 

To  aid  in  the  understanding  of  how  various  geometric  and  stiffness 
quantities  interact  to  effect  the  rod  motion,  several  more  specialized  cases 
will  next  be  examined. 

2.  Straight  Rod  With  Variable  Section  Properties. 

If  X  and  k  and  their  derivatives  are  set  equal  to  zero,  the  following 
equations  are  obtained: 
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where  superscript  prime  indicates  differentiation  of  vector  function  P  with 
respect  to  z. 

Proof : 


<\  -  P,  -  P  -  h  P'  +yr  P"  +  ...  , 

h2 

=  P2  -  P  =  -  h  P'  +  jy  P”  -  ...  , 
x  =  h3P'xP"  +  ...  ,  and 
|  CjXt  |  =  |  h3  (  jp'xP'  '|  +  ... 


Note  that  T  and  0  are  perpendicular  to  each  other  since  PWP'xP'')  =  0 

**If  z  is  the  arc  length  s  from  a  fixed  point  PQ  of  C  to  point  P,  0  is 
then  given  by 


0  - 


dP 

ds 


d2p 

ds2 


d2p 

ds2 


Proof : 


Since 
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Therefore  P'(s)  is  perpendicular  to  P''(s)  so  that 
|p’(s)  xP"(s)|  =  |P’  j  |P”  |  =  |P"(s)  |  . 


**The  plane  passing  through  P  and  normal  to  3  is  called  the  plane  of 
curvature  or  oscillating  plane  of  C  at  P . 
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Referring  to  Figure  A-2,  and  P~  are  points  on  C  either  side  of  point  P 
such  that  z  increases  by  an  amount  h  from  P  to  P^  and  decreases  by  the  same 
amount  from  P  to  Pj.  and  G>  are  position  vectors  of  P^  and  P2  with  respect 
to  P  respectively 


C 


Figure  A-2  . 

A  vector  binormal,  0  ,  of  C  at  P  is  defined  as 


C1  x  s 

8  =  lim  17  —  -- ~j*  .  (A-4) 

h+0  lCl  X  SI 

Again,  the  sense  of  0  depends  upon  the  choice  of  scalar  variable  z  used  to 
describe  the  position  of  P. 

If  P  is  the  position  vector  of  a  point  P  on  C  relative  to  a  fixed  point 
0,  the  vector  binormal  can  be  expressed  as 
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A  vector  tangent,  t  ,  to  C  at  point  P  is  defined  by 
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where  h  is  the  increase  in  z  from  point  P  to  P .  Note  that  the  sense  of  x  is 
not  unique  since  it  depends  on  the  choice  of  scalar  variable  z.  It  points  in 
the  direction  in  which  z  increases  along  the  curve. 

**The  vector  tangent  can  be  expressed  as : 
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where  P'  is  the  derivative  of  P  with  respect  to  z 


Proof : 
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Pj  =  P( z+h) 


By  Taylor's  series  for  vector  functions, 
P(  z+h)  =  P  +  P'h  +  P"h2/2!  +  ... 
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**If  the  scalar  variable  governing  the  position  of  P  on  C  is  the  scalar 
arc-length  displacement  of  P  relative  to  a  fixed  point  PQ  of  C  then 
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SPACE  CURVE  GEOMETRY 

This  appendix  defines  and  presents  relationships  among  various  vector  and 
scalar  quantities  used  to  describe  a  line  curved  in  three  dimensional  space. 
These  quantities  include  the  vector  tangent,  vector  principal  normal,  vector 
binormal,  the  radius  of  curvature,  and  the  curvature  and  the  torsion  of  a 
curve  at  a  point . 

Although  the  ensuing  material  may  be  found  in  many  sources,  this 
particular  exposition  is  abstracted  from  the  book  by  T.R.  Kane.*  Rigorous 
proofs,  which  have  been  omitted  from  this  presentation  in  favor  of  intuitive 
arguments,  can  be  found  in  that  book. 

1  .  Vector  Tangents  Of  A  Space  Qirve  (t  ) 

In  Figure  A-i ,  C  is  a  curve  in  space,  0  is  a  point  fixed  in  space  while  P 
and  are  points  on  C  whose  position  along  C  depends  on  a  variable  z. 

Vectors  P  and  ,  are  position  vectors  of  P  and  P^  respectively,  with  respect 
to  0  while  C  is  the  position  vector  of  with  respect  to  P. 


C 


Figure  A-l . 


T.R.  Kane ,  Analytical  Elements  of  Mechanics:  Vol.  2  Dynamics ,  Academic 
Press ,  New  'fork,  1981. 
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V.  CONCLUSIONS 


This  work  was  motivated  by  a  requirement  for  a  general  analytical  model  of 
a  gun  tube  which  could  be  used  to  examine  the  effects  on  tube  motion  of 
centroidal  axis  curvature  and  twist,  cross  section  eccentricity  and  variation 
in  cross  section  shape  along  the  axis. 

The  model  which  has  been  formulated  is  based  upon  first  principles  of 
mechanics  and  is  sufficiently  general  to  permit  study  of  the  above  effects  and 
also  of  motions  involving  large  displacements  and  rotations  of  the  tube  axis. 

The  model  can  provide  useful  information  in  its  own  right;  both  by 
illustrating  the  manner  in  which  the  governing  equations  are  coupled  when 
terras  representing  various  affects  are  added  or  deleted,  and  by  the  free  and 
forced  response  solutions  which  may  be  obtained  for  various  special  cases.  It 
can  also  be  used  to  provide  guidance  for  the  development  and  use  of  more 
approximate  analytical  techniques,  such  as  finite  element  models,  by 
establishing  the  potential  error  associated  with  various  approximations  to  the 
actual  structure.  Since  the  derived  equations  show  the  coupling  among  the 
displacement  and  rotation  variables  caused  by  axis  curvature  and  cross  section 
non-uniformity  effects,  they  can  aid  in  the  interpretation  of  physical 
measurements  of  gun  tube  motion. 
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Introduction  of  torsion  again  causes  coupling  among  all  of  the  equations 
nvun  for  a  rod  with  symmetric  cross  section.  It  is  also  noted  that  although 
the  equations  describing  motion  of  a  rod  with  constant  cross  sectional 
properties  would  have  somewhat  fewer  and  simpler  coefficients,  they  would 
retain  the  coupling  among  the  motions  described  by  the  four  variables. 
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It  is  seen  that  while  inclusion  of  effects  of  variable  curvature  somewhat 
complicates  the  coefficients  and  increases  their  number,  the  equation  with  I 
equal  to  zero  still  describes  two  separate  pairs  of  coupled  motion  involving 
vQ  and  $  together  and  uQ  and  wQ. 

5.  Hod  With  Constant  Curvature  And  Torsion  And  With  Variable  Cross 


Section 

The  final  special  case  of  Eq.  (11)  which  will  be  considered  is  for  a  rod 
with  constant  torsion  and  curvature  and  variable  but  axisymmetric  cross 
section  so  that  I  remains  zero.  The  displacement  equilibrium  equations  for 
this  case  are  given  below: 
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Eqs.  (13)  are  seen  to  be  fully  coupled  in  the  variables.  If  however,  I  is 

xy 

zero  as,  for  example,  for  a  circular  cross  section,  then  two  pair  of  equations 
result;  one  pair  involving  vQ  and  <J>  only  and  the  other  uQ  and  wQ. 

4.  Rod  Curved  In  A  Plane  With  Variable  Curvature  But  Constant  Section 
Properties, 

The  next  special  case  is  that  of  a  rod  for  which  is  zero  but  the 
curvature  completely  variable.  The  principal  normal  and  binomial  directions 
are  assumed  to  be  principal  directions  for  the  cross  section.  Eqs.  (11) 
become 
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Eqs.  (12)  are  seen  to  be  the  usual  equations  of  motion  of  shafts,  beams  and 
bars.  Equations  involving  uQ  and  vq  are  uncoupled  if  the  x  and  y  axes  are  the 
principal  axes  of  the  cross  section. 

3.  Rod  Carved  In  A  Plane  With  Constant  Curvature  And  Variable  Section 
Properties. 

In  this  case,  the  rod  is  initially  in  the  form  of  an  arc  of  a  circle  so 
that  kq  is  constant  and  is  zero.  The  governing  equations  thus  become 
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3.  The  Vector  Principal  Mor— 1  Of  A  8pace  Curve  (  v) 


The  unit  vector  given  by 


v  *  0  x  T  (A-7) 

is  called  the  vector  principal  normal  of  C  or  P.  Contrary  to  the  cases 
of  B  and  t  ,  the  sense  of  v  is  unique* 

An  alternative  expression  for  v  in  terras  of  derivatives  of  position 
vector  P  with  respect  to  coordinate  z  is 


(P'xP1 ')  x  P' 

|(P'  X  P")  X  P'l  * 


(A-8) 


In  terras  of  derivatives  of  P  with  respect  to  s.  v  is  given  by 


"(s 
P"(s)|  * 


(A-9) 


This  follows  from  the  fact  that  P''(s)  is  perpendicular  to  P'(s)  and 
P'(s)  is  a  unit  vector.  The  plane  passing  through  P  and  normal  to  v  is 
called  the  rectifying  plane  of  C  at  P, 


4.  The  Vector  Radiua  Of  Curvature  Of  A  Space  Curve  (p) 


Referring  to  Figure  A-3,  P  is  a  point  on  C  and  P,  is  the  point  of  C  with 
which  P  coincides  when  scalar  variable  t  is  increased  by  an  amount  h.  Cj  is 
the  position  vector  of  Pi  relative  to  P.  By  definition,  the  center  of 
curvature  of  C  at  P  is  the  point  approached,  as  P^  approaches  P  by  the  center 
Q,  of  a  circle  which  is  tangent  to  the  tangent  line  of  C  at  P  and  passes 
through  P2.  If  r  is  the  position  vector  of  Q  relative  to  P,  the  vector  radius 
of  curvature,  p  ,  is  defined  by 


p  -  Lim  r  .  (A-10) 

h+0 

**In  terms  of  the  derivatives  with  respect  to  coordinate  z  of  the  position 
vector,  P,  of  P  relative  to  a  point  0  fixed  in  space,  the  vector  radius  of 
curvature  is  given  by: 


(p» ) 
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Proof:  r(by  construction)  is  perpendicular  to  both  T  and  T  x  C. .  Hence, 
there  exists  some  scalar,  say  5  such  that 


r  *  St  x  (1  x  Cj)  , 


(A) 


On  the  other  hand,  as  Q  lies  in  the  perpendicular  bisector  of 

line  P  Pj  ,  r  may  be  regarded  as  the  sum  of  V2  Cj  and  a  vector  perpendicular 

to  both  Cj  and  txc^  .  Hence,  there  exists  some  scalar,  say  y  ,  such  that 


r 


fl 
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+  yCjX(x  X  Cj)  . 


(B) 


Dot  multiply  each  of  Eqs.  (A)  and  (B)  with  Cj  (to  eliminate  6),  equate  the 
resulting  expressions  for  r*^,  solve  for  5  and  substitute  into  Eq.  (A)  yields 


c^x  x  (t  X  Cj) 
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Next,  use  the  relationships 


r  *  TFT  :  e  '  h  p'  *  1,2/2  p 


and  substitute  Eq.  C  into  Eq.  A-10. 


5.  The  ladius  Of  Curvature. 


The  vector  radius  of  curvature  p  can  be  expressed  as 


P  *  P  v 


(A-12 ) 


where  p  is  an  intrinsically  positive  quantity  called  the  radius  of  curvature 
of  C  at  P  which  is  given  by 


t~D"T  • 


(A-13 ) 


Proof:  Multiply  the  numerator  and  denominator  of  the  expressions  for  P  of  Eq. 
(A-ll)  with  |p'|  ,  giving 


P'  _  (P'xP")  x  P' 

p  *  ‘  P'xP' '  I  Ip'xp"]  Ip' I  * 


By  Eq.  (A-8),  Eq.  (A-12)  is  an  identity  if  P  is  defined  by  Eq.  (A-13). 


**In  terras  of  derivatives  of  P  with  respect  to  8,  p  and  p  are  respectively 
given  by 


Proof:  |p* (a )  2  ■  1  and 


tP'(s)  x  P"(s)]2  -  |P'(s)  x  P* '  (a)  | 2  -  |p"(a)|2  . 


Furthermore, 


(P'(s)  x  P"(s))  x  P’(s)  -  (P')2P'*  -  P"*P'  P*  -  P"(s)  . 


6.  The  Serret-Frenet  Formulae 


The  derivatives  of  the  vectors  r, 


B  and  v  with  respect  to  s  are  given  bj 


dx 

ds 


v/p  , 


(A-16) 


d£ 

ds 


X  v  ,  and 


(A-17 ) 


dj» 

ds 


XU  -  x/p 


(A-18) 


where  p  is  the  radius  of  curvature  and  X,  called  the  torsion  of  the  curve  c 
point  P  is  given  by 


X  *  p2[P'(s)  •  P”(s)  x  P,"(s)]  . 


(A-19) 


Eqs.  (16),  (17)  and  (18)  are  called  the  Serret-Frenet  formulae. 
Proofs: 


at 
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T  -  P'(a) 

T'(a)  -  P"  -  P"  -  |P,,|  v  -  v/p  , 

8  -  P'(s)  x  P"(s)[(P,,)2]"1/2  ,  and 
B'(s)  “  P'xP,,'[(P,,)2f  1/2 

-  P'xP’'[(P")2r  3/2  P"«P,M  . 

Eliminate  the  first  and  second  derivatives  of  P 

P'  *  x  and  P"  *  v/p  . 

Then 


8'  =  p  t  x  (P'  1  '  -  v  vP'  1  1 ) 

=  P  t  x  [v  x  (P' ' 'xv)l 

■  P  T  •(P",XV)  V  -  p  T  •  V  (P"'XV)  . 

Since  , 

V*  T  ■  0, 

8'  *  pt  •  (P' "  x  v)  w 
-  -  X  v  , 

where 

X  *  -  p  T  •  (P"  'XV) 

»  P2  [P*  •  P"xP'  "  J  . 


(A-9) 


by  using 
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-  -  ,\  ,t 


Finally , 

v  ■  0  x  x  ,  and 
v'  -  0'  X  X  +  0  X  X1 

*-\VxX+0X  V/p 

But  • 


v  X  x  *  (0xx)  x  x  =  -  0  and 
0xV  -  0x(0xx)  =  -  t  . 


Hence  , 


w'  -  A0  -  x/p  . 


**The  tors  ion, A,  expressed  in  terms  of  derivatives  of  P  with  respect  to  a 
variable  z  other  than  8  is  given  by: 


A  »  p2  [P'(z)  •  P”(z)  x  P'"(z)]|p,(z) 


(A-20) 


Proof: 


P'(s) 


> 


P"(s) 


p"(«)  (If-)2  *  »’<*> 


ds 


and 
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»  P-  »  * laT*1  ,  a  *  fc  P  i  ^  "  _*!, 


+  p' 


P’"<2)(£)3  +  3  P"(z>  ^4- 


Then  , 


P*(s)  •  P”(s)xP'”(s)  =  P'(z)  •  P,,(z)xP,,,(z)(—)6  , 


(&  P'(z)2  =  P’(s)2  -  1  . 


Therefore, 


r  dz\6 


[P’(z)2]  "3 


P'(z)  ' 


7.  Example ;  The  use  of  the  Serret-Frenet  formulas  is  illustrated  by  use  of 
the  following  example . 

A  straight  line  AC  is  drawn  on  a  rectangular  sheet  of  paper  ABCD  having 

the  dimensions  shown  in  Figure  A-4  .  The  paper  is  then  folded  to  form  a  right 

cylinder  of  radius  a/ir,  the  line  AC  thereby  being  transformed  into  the 
circular  helix  H  shown. 

a.  Determine  the  cosine  of  the  angle  ij>  between  a  unit  vector  k  parallel 
to  the  axis  of  the  cylinder  and  the  vector  binormal  of  H  at  point  P. 

Solution:  Let  P*  be  a  point  on  the  cylinder  base  (Line  AB)  such 
that  PP'  is  parallel  to  k.  Let  z  be  the  distance  between  P  and  P',  0  the 

intersection  of  the  axis  of  the  cylinder  with  its  base,  n  a  unit  vector 

parallel  to  line  OP',  9  the  radian  measure  of  angle  AOP',  and  P  the  position 
vector  of  P  relative  to  0. 


‘.a  v  V  V  *  ■  V  *  -  *  •  *  *  V  V  *  •  ‘  *  •'  v  •  •  .  •>  ■ 


P'  ■  a/n  k  x  n  —  +  k  , 


8  **  n  —  (Since  arc  AP'  ■ 


z  —  ,  and 

D 


d8  it 

dz  "  b  * 


Therefore . 


P'(z)  =  -g  kxn  +  k  , 


'"<*>  te(kx.)  £ 


~  ir  n  (Since  kx(kxn)  **  -  n)  . 
b 


P'xP"  =»  ■—  k  +  n  x  k)  and 


|p'xP'  ’I  »  it  ~2  ( 1  +  1/2 


2  2 ' 

b  o 


Thus  , 


0  *  (-  k  +  n  x  k)(l  +  — )‘ 


Finally 


cos  Ip  -  k*  0-  ( 1  +  ^2  . 

& 


b.  Determine  The  Radius  Of  Curvature  At  P  Of  The  Helix  H. 


c.  Determine  the  torsion  \  of  the  curve  H. 

A  -  p2[P,»P,,xP,,')[P']’,S  , 

P'  ■  —  k  x  n  +  k  , 
b  ’ 


pi  i  • 


dm 

dz 


2 

a  % 

x  —  kxm  , 
bJ 
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•  Cm 


Then  , 


M  -  (i  +4) 1/2 » and 


P'*(P"xP"  ’)  -  l  a-r-  . 

b5 


2^.2 
a  +b 


zular  Velocity  Of  Space  Curve  Coordinates  In  A  Fixed  Reference  Ftase. 


Let  P  be  a  point  moving  along  a  space  curve  c  with  speed  —  where  8  is 
the  distance  to  P  along  c  from  some  point  of  c .  If  R  is  a  reference  frame  in 
which  the  tangent  (t),  principal  normal  (v)  and  binormal  (0)  unit  vectors  to  c 
at  P  are  fixed  then  the  angular  velocity  of  R  in  the  reference  frame  R'  in 
which  c  is  fixed  is  given  by 


R'uR»ai  +  i/p  u)g. 


(A-21) 


Proof:  Given  two  reference  frames  R  and  R'  and  any  two  non-parallel  vectors  a 
and  b  fixed  in  R.  The  angular  velocity  of  R  in  R'  is  defined  by  Kane  as 


R'  R 
to 


da  db 

dt  X  dt 
da  . 

3?  •  b 


(A-22) 


If  a  and  b  are  taken  as  the  vector  principal  normal,  »  ,  and  the  vector 
binormal,  P,  respectively  at  P  then  Eq.  (22)  may  be  written  as 


dw  dg 

ds  X  ds  ds 


(A-23) 


Next,  the  Serret-Frenet  formulas  give 


dv 

ds 


XB  -  1/p  T  and 


d$ 

ds 


X  v  , 


Substitution  of  the  above  into  Eq.  (A~22)  gives  the  desired  result. 


I 

i 
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APPENDIX  B 


APPENDIX  B.  COEFFICIENTS  APPEARING  IN  EQUATIONS  ( IV- I 1 ) 
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D11  ■  <GV 
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